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... Abstract. The nonlinear Schrodinger-Helmholtz (SH) equation in A'^ space dimensions with 2(7 non- 

' linear power was proposed as a regularization of the classical nonlinear Sehrodinger (NLS) equation. 

C " 3 ^ It was shown that the SH equation has a larger regime (l<CT<-^)of global existence and unique- 

■ ness of solutions compared to that of the classical NLS (0 < (t < -^). In the limiting case where the 
(N ■ Schrodinger-Helmholtz equation is viewed as a perturbed system of the classical NLS equation, we apply 

^ modulation theory to the classical critical case (cr = 1, A'^ = 2) and show that the regularization prevents 

, the formation of singularities of the NLS equation. Our theoretical results are supported by numerical 

■ simulations. 
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1. Introduction 
The nonlinear Schrodinger-Helmholtz (SH) system is given as below 



(N 

u~a^Au=\v\''+\ (1) 

■ v(0)^vo. 

■ 

where cr > 1 and, for simplicity, a > 0. This system has been proposed [5] as a regularization for the 
classical nonlinear Sehrodinger (NLS) equation: 

00 ! iut + Au+\u\^''u = 0, t>0,xeR^, 

'. u{0) = UQ. (2) 

• ^ . In [5] we showed global existence of solution of the Cauchy problem ([1]) for 1 < cr < 3 when N = 1 and 
' 1 < cr < when > 1. It is well-known that the classical NLS has global solution for < cr < in any 

■ dimension > 1 and there is finite time blow up in the critical case a ~ ;^(see e.g., [5], [12], [II], [T5] . 
[TO] , [5D] and references therein) . So we regard the SH system ([T]) as a regularization system for the NLS 
([2]) since the former system has larger regime of global existence for the parameter cr, which contains the 
values for which the NLS ([2]) blows up. Note that ([T]) is a Hamiltonian system with the corresponding 
Hamiltonian 

H{v) = / I V^;(x, t)p - ^ ' ^' ^ ^' dx 

jR" V cr + 1 y 

and can be obtained formally by the variational principle 

.dv sniv) 



i- 



dt Sv* ' 

where v* denotes the complex conjugate of v. Let us rewrite system ((1} as 

ivt + Av + Iwp'^u + a'^{Au)\v\''-'^v = 0, (3) 

where u = (/ — Q!^A)~^(|w|'^+^). Observe that when the parameter a goes to zero, one can regard system 
(H]) (or [3]) as a formal perturbation of the classical NLS. There has been a lot of work on perturbed 
NLS in the critical case a = jj { see, e.g., [S], [3], [TU], and references therein). In this paper, we will 
apply modulation theory (see e.g., [TT], [T^], [T7] and [H] for references about modulation theory) to 
the classical critical case cr = l,Af = 2 in order to shed more light on the nature of the effect of the 
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regularization in preventing the blow up. In this case the classical NLS blows up for certain initial data, 
however, the SH system has global solution with the regularization parameter a > 0. Indeed, modulation 
theory tries to explain the role of the regularization in preventing the formation of singularity near the 
critical values of the initial data which blow up in the classical case. Intuitively speaking, the basic 
idea behind modulation theory is that the energy near singularity is equal to the power of the Townes 
soliton [see Eq. ^ below] , and the profiles of the solutions are asymptotic to some rescaled profiles of the 
Townes soliton. With modulation theory, one can reduce the perturbed system ([3]) into a simpler system 
of ordinary differential equations that do not depend on the spatial variables, and they are supposed to 
be easier to analyze both analytically and numerically. 

In this paper, we will study how the parameter a prevents the singularity formation in the 2D critical 
NLS. The work here will follow the study we initiated in [5]. In particular, some of the statements have 
already been mentioned there. For the sake of completeness, we will restate some important theorems 
and propositions. 



2. Modulation theory 

First we review some main results on modulation theory for the unperturbed critical NLS following 
[l2] . As stated in [12] , most of the results presented in this section are formal and have not been made 
rigorous at present. We emphasize here that we consider the case cr = 1, TV = 2 in order to see how the 
regularization prevents singularity formation. 

In the case of self-focusing the amount of power which goes into the singularity is equal to the critical 
power Nc = \\R\\2 ~ R^{r)r dr where R, the Townes soliton, is the solution of the following equation 
with minimal L^— norm, which is positive and radially symmetric 

AR-R + R^ = 0, i?'(0)=0, limi?(r) = 0, A = — + - — • (4) 

i-^+oo dr-' r or 

Close to the stage of self-focusing, the solution of ([2]) separates into two components as it propagates, 

v^v,+ «back (5) 

where Vs is the high intensity inner core of the beam which self-focuses toward its center axis and ^^j^^ck 
is the low intensity outer part which propagates forward following the usual linear propagation mode. 
Close enough to the singularity, Vs has only small excess power above the critical one and approches the 
radially symmetric asymptotic profile: 

Vs{x,t)^^ViT,Ocxp(^i{r+^^)y argViT,0)^0, (6) 

where L{t) is a yet undetermined function that is used to rescale the variables ^ = (■^1,^2), x = {xi,X2) 
with 

X dr 1 
L' dt ^ L^' 

Then the reduced system for unperturbed critical NLS is 

where is the scaling factor and /? is proportional to the excess of the power near singularity: 

/3 = M{N — Nc) for constant M = j Jq R^{p)p^ dp ~ 0.55. We emphasize here that in the case of 
self-focusing of the original system, both the scaling factor L and the excess energy /3 will approach zero. 

Next we review some results on modulation theory for the perturbed critical NSL [12]. For a general 
perturbed critical NLS of the form 

ivt + Av+\v\'^v + eF{v,vt,Vv,- ■ ■) ^0, |e| < 1, (8) 

where F is an even function in x, modulation theory is valid when the following three conditions hold. 
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Condition 1. The focusing part of the sohition is cfose to the asymptotic profile 

Vsit,x) ^ -i-^F(T,Oexp[*T(t)+z^^], (9) 

where 

and V = R. + 0{(3, e), /? = —L^Lu and R is the Townes soHton given in 
Condition 2. The power is close to critical 

I 2^ J \v,{t,Xi,X2)\^ dxidx2- N^\<.1, (10) 

or, cquivalcntly, 

1/3(01 «1, (11) 
where = Jjja R{xi,X2)^ dxi dx2 = \\R\\2 is the threshold energy of blowup. 

Condition 3. The perturbation eF is small in comparison to the other terms, i.e., 

\eF\ < |At;|, |eF| < \v\\ (12) 

The following proposition is given in |12| . 

Proposition 1. // conditions 1-3 hold, self-focusing in the perturbed critical NLS (0) is given to leading 
order by the reduced system 

/3t + ^ = ^(/i)t-^/2, Lu = ~^. (13) 
The auxiliary functions fi , /2 are given by 



flit) = 2L{t)Re 



f2{t) = Im 

where 



^ I F{^r) eM-iS)[R{p) + pVi?(p)]rfxidx2 
I ip*RF{'4;u)dxidx2 



(14) 
(15) 



i'R = -^■R(p)exp(iS'), 



and R is the Townes soliton given in p = j-, S = r(t) + -^^, ^ ~ M = j R{pY dp ~ 0.55. 

Furthermore, if F is a conservative perturbation, i.e., 

Im / v*F{v)dxidx2 = 0, (16) 

then f2=0. Since /3 ^ 1, e <^ P, taking the leading order by neglecting the exponential term in the 
first equation of 113\) . we further reduce the system I113\} into the following system 

-L^Lu=f3^Po + ^fi. /3o=/3(0)-^/i(0), (17) 
where is independent oft. 

In general, at the onset of self-focusing only condition 3 holds. Therefore, if the power is above Nc the 
solution will initially self-focus as in the unperturbed critical NLS. As a result, near the time of blowup 
in the absence of the perturbation, conditions 1 — 2 will also be satisfied. 

It is worth pointing out that, as studied in |12| . various conservative perturbations of the critical NLS 
equation, for instance, self- focusing in fiber arrays (see [T], [3], [1], (TB], [H] and references therein) 
and small dispersive fifth-power nonlinear perturbation to the classical NLS |18j . have a generic form 

C 

/i — — , C = constant, (18) 
which results in a canonical focusing-defocusing oscillation. 
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Next, we shall derive the reduced equations (fT7|) that corresponds to the nonlinear Schrodinger- 
Helmholtz regularization system in the critical case cr = 1, N = 2 (the reason for this restriction is that 
it allows us to compare with numerical simulations). In this case, Eq. ([3]) reads 

ivt + Av + \v\'^v + a'^vAu = . (19) 

Comparing Eq. (fT9|) with ([8|) we have e = and 

F{v)^vAu, li- a^Aw |t>p. (20) 

We shall assume that the system ([T^ satisfies all three conditions, ([5]), (UHl) and p^ . Since 

^i?(a^) = ;^i?(|)cxp(*^), (21) 

we have 

F{i;R){x) = i;Rix)AuR, (22) 

where ur satisfies 

UR{x)^a^AuR{x) ^ \^jr{x)\^ = |ii?(|)exp(z^)|2 = -^|i?(|)P . (23) 
For a given function 5, solution to the equation 

is given by 



(I -a^A)u{x) ^ gi-), a; G (24) 



uix) = iB..g)i^), 



where 5° is the modified Bessel potential or the Green function corresponding to the Helmholtz operator 
(see e.g., [7] for reference on Bcssel potential) 



1 1"^ e'^e 

Bfi^) = ^ Zm-. ds. (25) 



2a2 sAf/2 
If we let g{-) ~ ■^R{-)^ then we can write the solution to Eq. as 



ob(x) = (i3^,_R2)( )= (Bj,r2)(_) (26) 



1 1 f I e-^e 
/r2 I JO 



y ; ds R\y)dy^dy2, (27) 



where y = {yi,y2)- Substituting the above into ([^ and using we get 

F(V'fl)(x) = Vi?(a^)Aui?(x) = ^(iii^(a;)-|Vi?(x)|2)^fl(x) (28) 



ur{x)~^R\^) 



(29) 



Now we can calculate the term /i 

A = -i?e / F(?^i^(x))exp(-^5)(i?(p)+pi^p)da;lda;2 

L f 1 1 

= -/ ^^(/o)[ui?(2;)- Ti^^(p)][^(p)+P^p]da;irfa;2 

= ^-7 / Ufl(a;)i?(/9)[i?(/9) + /9i?p]da;idx2 -^t^ / i?^(p)[i?(p) + pi?p]da;ida;2 

= J1-J2. 

It is easy to see that the second integral J2 is a constant that does not depend on L. Indeed, applying 



change of variables: ^1 — = 77 ' since P = j = l ~ V^i + ^2; S^t 

J2 = ^ / R\p){R{p) + pRp{p)) diid£.2 = CO, (30) 
where Co is a constant that does not depend on L. 
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Next, let us look at the first integral Ji. Plugging the result of (|27|) into Ji, we get 

Ji^-^ ur{x)R{p){R{p) + pRp)dxidx2 (31) 
7ra^ Jr2 

ds I i?"(y) dyidya R{p){R{p) + pRp)dxidx2 (32) 



ds i?2(y) dyidy2 R{p){R{p) + pRp)dxidx2. (33) 




Change of variables: ^ = (€1,^2) = (7^1 7^)1 then we will have 



L ' L 



J^=o 4/ / / dsR\y)dyidy2R{p)(R{p)+pRp)d^id^2- (34) 

So /i can be written as 



fi = 7y 1 / / / ds R\y) dy,dy2 R{p){R{p) + pRp)d^i - cq, (35) 

27r Jr2 Jr2 Jo s 

where cq is a constant given in (|30p . 

Plugging into (fT7|) we have the reduced system for the Schrodinger-Helmholtz system ([T|) 

2 2 
- L^Lu = /3o + ^/i, /?o = m ~ (36) 

where /i is given in ((35|) . 

Now one needs to study the ordinary differential equation (|36p with /i given in ([35)l . The explicit 
form of the function /i is much more complicated when compared to the generic form of p8|) due to the 
nonlocal nature of our perturbation term . The idea is to show that for L small this additional term 
on the right-hand side of the first equation of ([36]) will prevent the singularity formation, i.e., prevents 
L from tending to zero as time evolves. In the next section, we will investigate the ODE system pop by 
approximating the function fi . 

3. A SIMPLIFIED Reduced system 

Now, without solving Ui?(x) explicitly as we did above, we will try to approximate ur{x) by asymptotic 
expansion in terms of ^ , for small values of ^ , and further approximate the function /i . 

3.1. First order approximation. Recall that from ([T]) we have (when (7 = 1, iV = 2) 

u — a^Au = \v\'^ 

or we can write u{x) = (/— a^A)^^\v{x)\'^ . When a is very small, we can formally write u{x) in the first 
leading order term: u{x) = \v{x)\'^ + 0{a'^). Now for iPh{x) = j^R{j^) exp {iS), we can similarly write 
ur{x) as 

Ur{x) = \i;R{x)\' + 0{{^f) 

= I2^'(f)+0((f)') (37) 

so we have 

Fii^R)ix) = {A,,UR{x))i;R{x) 

^ {A,i\^R{x)\^))i;R{x) 

^ -l(A.i^2(|))ii^(|)exp(^5) (38) 
Substituting this into the equation for /i p^ . we get 

/i = -Re [ F{i;Rix))expi^iS){Rip)+pRp)dxidx2 

^ -/ ±{A^R^^))iR{^){R{p) + pRp)dx,dx2 (39) 
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Next, we make change of variables: 



then for p = |^|, and by the chain rule, we get 



L 



2 h (40) 



where 



/ {A^R^)R{R + pRp)d^id£,2 

/"OO 

= 2 / \{R^)pfpdp > 0, (41) 
Ja 

where the detail of the calculation of the above integral is presented in claim 1 of the appendix. 

Plugging into (jSG]) . we have the leading order of the reduced system, which turns out to be of the 
generic form 

2 

- L^Lu = /3o + ^/i (42) 

with /i ^ and Pq = /3(0) + ^^x^^o) < 1 since /3(0) < 1, f < 1. Fibich and Papanicolaou [H] 

showed that there is no singularity in finite time with this perturbation. In fact, substituting /i 
into the above equation, we get 

Write y = L^, then y satisfies the following equation: 

{y,f=4(3o-^-+4Doy, (44) 
M y 

where Dq is a constant satisfying Dq ~ (0) — T^fs) + "aa/ l^o) ■ 

Since —jg^ > 0, y can not go to zero in the above equation, i.e., L can not go to zero, which explains the 
prevention of the singularity formation, at this leading order in the expansion. 

3.2. Second order approximation. In the previous subsection we use the asymptotic expansion by 
taking the first leading term as approximation for the solution u of the Helmholtz equation u — a^Au = 
One might naturally ask whether we will get better approximation and still have the no blow up 
structure if we approximate the solution of the Helmholtz equation u{x) by taking one more leading 
term. To investigate this, we proceed similarly as before: we can formally write u(x) in the first two 
leading terms: u{x) ~ + A\v{x)\'^ + 0(q;^). so we have 

URix) = |Vi?(x)p+a2A,|^«(a;)p+0((^)4) 

= -li?2(|) + a2A.(-li?2(|)) + 0((^)4) 

- ^,R\l) + A^R^l) (45) 



and then 



F{,i1^r){x) = {A,ur{x))^r{x) 



A.(-^i?^(|) + (|)2A.i?2(|))) cxp {^S) 

-lA.i?2(|) + qfAlR\l)\ li?(|) exp {^S) 
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Substituting this into the equation of /i ()14|) . we get 

/i = -Re [ F{ijR{x))exp{~iS){Rip)+pRp)dxidx2 

- {]^A,R^i^) + i^)^AlR^i^))^Ri^){R{p) + pR,)dx,dx, 
~ - / ±{A.,R\j))R(j){R{p) + pRp)dxidx2 

+ - f {jfAlR^{j)R{j){R{p) + pR,)dx,dx2 

- /l + /2 

Here the first integral Ii is equal to — ^ with Ci in (|^T|) as we calculated in the first order expansion 
case. For the second integral I2, after change of variables C = = (^11^2) and p = |^|, it gives 



h = 



/ {Al{R^))R{R + pRp)di,di2 



where 



(46) 



C2 = - I iAW))R{R + pRp)d^id^2 

= f iA^R' fd^,d^2 > 0, (47) 
where the detailed calculation of the above integral is presented in claim 2 of the appendix. 
Therefore, we obtain 

h--Jl + ^ (48) 

with a correction next order term to the generic form (|18|) . It yields the reduced equation of the leading 
order 

-L'Lu = f3o + ^fi. (49) 
Substituting /i into the above equation ([49|) . we get 

Then y = Li^ satisfies the following equation: 

/ n2 a^C*! 1 2a4C2 1 ^ ,,,, 

(y*)' = 4/3o - -r^- + + E^y, (51) 



where -Bo is a contant satisfying £'0 = 4i^(0) - + ^x?(o) ~ W 

Let us rewrite the right-hand side of the above equation (|?T|) by substituting y = and i?o into the 
equation 



When L approaches zero, is very small and in the regime of <^ 1, which is required assumption 
for the expansion, the right-hand side will remain positive when Lt{0) is large. In other words, contrary 
to the equation of first order expansion (|44p . the second order expansion equation (j50p might blow up 
with certain large Lt{0), we will see this is also verified in the numerical computation in the next section. 
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3.3. Third order expansion. With finite time blow up in the second order expansion, one miglit try 
to include a higher order term in the expansion of the solution of the Hclmholtz equation: 

u{x) = \v{x)\^ + a^A\v{x)\^ + a^A^\v{x)\^ + Oia"^). (52) 

As a result we have 

Unix) = + (f )^A.i?^(|) + qfAlR^l) + Oii^n 

and then 

F{iPr){x) = (A,Ufl(a;))Vfl(a;) 

A,.i?^(|) + |J A^i?2(|) + ^ A^i?2(|)) exp {^S). 

Now /i satisfies the following expression 

/i = -Re I F{i;Rix))exp{-iS){Rip)+ pRp)dxidx2 
- h + h + h, 

where Ii, I2 are the same as in the first and second order expansion ([^0]) and and we have 



h = 



- I ^[AlR\j))\R[j){R{p) + pR,)dx,dx2 
[ iAlR^)RiR + pRp)d^,d^2 



la* 



(53) 



L6 ' 

where the constant 

C3 = - / {\/AR')'d(,d^2 > 0, e = 7 = (a, 6), (54) 
where the detailed calculation is presented in claim 3 of the appendix. 



Substituting /i into the reduced equation of L (|42|) . we get 



So y = L satisfies the following equation 

where Fq = ALj (0) - -p^ + - 1 ^^^^ + . In this equation, since C3 is positive, we 

can see again that y can not approach zero, i.e., L can not go to zero, which prevents singularity formation. 



4. Numerical Results 

In this section we will show some numerical results of the evolution of L{t) in three expansion cases: 
(|43|) . (|50p and ([55]) in order to study the prevention of blow up of the Schrodinger-Helmholtz system. 
We will first consider the first order expansion case. 

First order expansion 

Let us look at the equation (|43|) . After some algebraic calculation and calculus integration, we come up 
with equation (|44)) : 

. .9 o?Ci 1 

ytf =4/3o--rTi-+47^oy, 
M y 

where /3o = /3(0) + and D, = _ ^ + l^^^" 

^From above we know that y can not approach zero, equivalently, L can not go to zero, which prevents 
singularity formation. 
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Furthermore, Fibich and Papanicolaou [T2] derived the generic equation 

i^^Cil Ho -4go 1 , 

h 4 y = 

M y My 



where 



{ytf = 4/3o - ^^T^- + 4^2/ = -^^-{vm - y)(2/ - ym), 



= ^21/07^17 = + 



J?o = H(0) 



2M ^/32 + a^CiHo/AP + /3o 4A//3o^ 



4 £4(0)- 



^From above we see that when Ho > 0, and Lt{0) > 0, i is monotonically defocusing to infinity; when 
Ho > and Lt{0) < 0, self-focusing is arrested when L = Lm ~ (j/m)^^^ > 0, after which L is monotoni- 
cally defocusing to infinity; when Ho < 0, then L goes through periodic oscillation between Lm = {ymy^"^ 
and Lai = {yAiY^^ (see Figure 1). 



Second order expansion 

Now wc look at the equation (|50|) of next order expansion: 



fC2 \ (^Y_ (Ci_ 



(f)'+/3o. (57) 



C 1 

In this equation, when (3o > c^gj/j the right-hand side is definitely positive, then for any initial data 
L(0) with Lt{0) < 0, the solution L will monotonically decrease and approach zero in finite time. In 
other words, when the initial excess energy is larger than certain amount (p^g^v/) ^^'^ ^ is initially 
focusing, L will focus and blow up in finite time (see Figure 2). However, this is not valid to begin with 
applying the modulation theory. Recall that for us to apply the modulation theory to a perturbed crit- 
ical NLS, we require three conditions to hold. One of the conditions is to require \(3{t)\ <^ 1 (see[TUlor[TT|). 



So we will consider the case when < /3o < c^g^i then we have r^^^ = y '"^ 2^02 ' ''high ~ 

'^^^2^C2^^ ' ^i^ere K ~ ^8l7 ^ Po- In this case, when initially > rj^jgj^ and £t(0) < 0, the 
right-hand side of the equation (j57[) will remain positive, so L will monotonically decrease to zero, which 



c 1 

is similar to the case of Po > chsM' ^^^'^^ again, this is not valid here for the discussion since the 
asymptotic expansion (|45p is valid under the assumption that ^ <C 1, so we need only to consider the 
situation of small, in this case, < ^^^[g]^- 



c 1 

Finally, we consider only the case < /3o < c^-gjj and < < r-^igii- 

First we focus on < < riow In this case, L might defocus to infinity, oscillate between two 
values or even blow up in finite time depending on different initial condition of Lt{0) for given /3o,Q! 
and L{Q). In Figure 3, we take the parameters a = 0.01, = riow/^. and /3o = 0.01. When initially 
Lt(0) > —44.9999, L will eventually defocus to infinity if Ho > (3a) and L will oscillate between two 
values if i/o < (3b); when initially Lt{Q) < —49.9999, L will approach zero in finite time, i.e., we 
observe singularity in finite time (3c). This numerical result is expected from analysis at the end of 
subsection p.2p . 

Similarly, for riow < -^o) ^ fhigh, we will observe different behaviors - defocusing, oscillation or fo- 
cusing depending on Lt(0) and Ho- For instance, when = , we have the threshold value 
X( = —39.9999, i.e., when Lt(0) > —39.9999, L will eventually defocus to infinity if Ho > and oscillate 
between two values if Ho < 0; when Lt{0) < —39.9999, L will eventually decrease to zero. 
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Third order expansion 

Lastly we study the equation ()55|) of higher order expansion: 



By defining y = and integrating the equation, we get equation ([5; 



My 6 M 2M 



wnere i^o - "i^t l^j - ^^(o) + TT TJ{o) ~ 3 UJo) + ^aT ■ 

In the above equation (|56p. y can not approach zero since the leading order term on the right-hand side is 
of negative sign when y goes to zero, equivalently, L can not approach zero. The nature of this equation 
is the same as that of equation (|44|) . and we see the same pattern in the numerical result (see Figure 4). 



5. CONCLUSION 

^From the analysis and numerical computation, we see that the regularization of the classical NLS 
effectively prevents singularity formation with positive parameter a > 0. By asymptotically expanding 
the solution of the Helmholtz equation to approximate the reduced system of the modulation theory, we 
observe strong no blow-up pattern in both first order and third order expansion. In the valid regime of the 
expansion and modulation theory, we also observe no blow-up pattern in the second order expansion with 
further restriction on certain condition: Lt{0) > L^, threshold initial value of Lt{0). This phenomenon is 
expected for even higher order expansion, say fourth order expansion approximation. One of the reasons 
is that the Laplace operator is not bounded, which causes instability for the expansion of the solution of 
the Helmholtz equation. 
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First order,alpha=0.01 ,beta=0.01 .Hp>O,L|(O)>0 




5 10 15 20 25 30 35 40 



la 



First order,aipha=0.01 ,beta=0.01 .Hp>O,L|(O)<0 




5 10 15 20 25 30 35 40 



16 

First order,alpha=0.01 ,beta=0.01 .H^<0 




Ic 

Figure 1. L evolves in time in first order expansion, la monotonic defocusing, Hq > 
0, it(0) > 0. lb first focusing then defocusing, Hq > 0, Lt{0) < 0. Ic oscillation, Hq < 0. 
For all cases, a = 0.01, Po = 0.01. la, 16, = 0.001. In Ic, = 1/8. 
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For completeness, we present in this section the detail of the calculation of the integrals (|40|). ([46|) and 
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Second order blow up 




Figure 2. Flow of L in second order expansion. L decreases to zero in finite tim, 
a = 0.01, /?o » 1, L{Q) = 0.1, Lt(0) = -2 



Claim 1: The integral Ii in (|40|) for the first order expansion can be simplified as 
1 1 f 2 

^ JR2 -tJ Jo 

Proo/. 

/i = --1/ {A^R\p))R{p){R{p)+pRp)d£,,d£,2 

= --^ / iAR^)RU^id^2 + -^ f iAR^)RRppdCidC2 

= hi + lu- 
lu the rest of this section, all the integrands and integrals are of variables ^ or p (i.e., they are the scaled 
variables) unless it is stated otherwise. 

Now for III, we integrate by parts once and change the variables by = pcos6, ^2 = psin6, we obtain 

1 1 



hi = - 



7rL2 



(Vi?2)2 d^2 



1 1 f°° 

= —-^2n / [{R')prpdp 

which gives us exactly the constant Ci as in (|¥T|) . 

For /12, we show now it is identically zero. 
Using polar coordinates, we get 

11/"°° 1 
/12 = -732^/ {{R^)pp + -{R^)p)RRpP^dp 



2 9 /"^ 

— / {R^)ppRRp P^dp+^ I {R^)pRRpP dp 
^j^ {R^)pp[R^)^p^dp+— j {R^)p{R^)ppdp 



L2 

X + Y. 
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Second order,alpha=0.01 ,beta=0.01 .alpha/L(0)=r /2, H >0 




5 10 15 



25 30 35 40 



Second order,alpha=0.01 ,beta=0.01 .alpha/L(0)=r 12, H <0 




25 30 35 40 



Second order,alpha=0.01,beta=0.01,alpha/L(0)=r 12 



5 6 7 



3a 



36 



3c 

Figure 3. Flow of L in second order expansion. 3a L dcfocuses to infinity when Hq > 0. 
35 L oscillates between two values when Hq < 0. 3c L blows up when |Lt(0)| is large. 

For all cases, we use a = 0.01, /?o = 0.01 << , a / L{0) ~ riowf^- 



For the first integral X, after rewriting and integration by parts once, we get 

1 r°° 1 



-Y. 
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Third order,alpha=0.01, beta =0.01, H||>01|(0)>0 




5 10 15 20 25 30 35 40 



4a 



Third order,alpha=0.01, beta =0.01, H||>01|(0)<0 




5 10 15 20 25 30 35 40 



46 

Third order,alpha=0.01 ,beta=0.01 ,Hq<0 
0.55 1 , , , , , , , 1 




■ 5 10 15 20 25 30 35 40 

4c 

Figure 4. Flow of L in third order expansion. 4a monotonic defocusing, Hq > 
0, it(0) > 0. 4& first focusing then defocusing, Hq > 0, Lt{0) < 0. 4c oscillation between 
two values, Hq < 0. For all cases, a = 0.01, /3o = 0.01 In 4a and 46, = 0.001. In 4c, 

We conclude that I12 = 0, which yields the result of and (|¥T|) . 

Claim 2 The integral /2 (|46p in the next order expansion can be simplified as 
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Proof. 

1 ^ /" 

h = / {A\R^))R{R + pR,)dC,d^2 



-^J iA\R^))R^d^,d^, + -^J iA'iR'))RipR,)d^,d^2 



— I2I + l22- 

After integration by parts twice, we get 

„2 



hi^-^J [A{R^)fdi^d£,2. 
Next we will show that I22 = 5-^21, which yields (|47|) . 

Recall that pRp = ^ ■ Vi?, ^ = (^1, ^2), = £,1 + ^o we have 

/22 = / (A2(i?2))(i?(^.Vi?))de 

^ -fj (A^(i?^))(e-vii?^)de 

= ifj/ [A(i?2)Pde + -fJ / A(i?^)(e-V(Aii?^))de 
= P + Q. 

For the second integral in the last line, we rewrite it then integrate by parts and obtain 

= 1^/ A(i?2)(e-v(Ali?2))de 

TT Jr2 2 

= / ^[A(«')l'<if 

2 

So we get /22 = -P + Q = = 5^21, so we conclude that I2 = ^hi ^ /r2 [A(i?^)]^ d^, which yields 
exactly (gH) and (gT]). " □ 



Claim 3 The integral (|53p in the calculation of higher order expansion can be simplified as 
^3 = -7? / (A|i?2)i?(i? + pRp) dC,d^2 = --^ / {VARY 



Proof. 



h = [ {A'{R'))R{R + pRp)d^ 

= f iA'iR'))R'd( + -^ [ iA'iR'))RipR,)d£ 



~ /31 + /32- 

After integration by parts three times, we get 



hi---^ I {VA{R^)fdi. 



16 Y.CAO, Z. H. MUSSLIMANI, AND E.S. TITI 

Next, we will show that I32 = hi = -^f^ /^^ (VA(i?2))2 ^j^, 
I32 



f iA'{R'mpR,)d^ 
I {A\R^))R{i-VR)di 
{A\R^m-w\R')di 



la* 

„4 



/ (A2(i?2))A(C.V-E2)de 



If^ / (A2(i?^)) f A(i?2) H 

(A2(i?2))A(i?2) 



e-v(A-i?2) 



A + B 



la4 
7rL6 



-(A-'{R'))i^-V{AR'))dC 



For ^, after integration by parts once, we obtain 

For i?, we define = A{R'^), a scalar function. Then we rewrite and calculate the term B as follows 



B = 



1 

7;-T7r27r 
27r L6 



a 



i^PP + -(t>p){p (t>p)pdp 



je I (t)pp(f>pP'^ + i(f>pfpdp 



0. 



So we have now I3 ~ I31 



^32 = 2/31 = -|fF /R2(VA(i?2))2 which concludes our resuh of ^ and 



□ 
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